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Abstract 

Using a Ginsburg-Landau free energy functional, we study the H — T phase 

diagram of a weak magnetic superconductor, where the magnetization from 

the magnetic component is marginal in supporting a spontaneous vortex phase 

in absence of external magnetic field. In particular, the competition between 

the spiral state and spontaneous vortex phase is analysed. Our theory is 

applied to understand the magnetic properties of EViV^-E^C- 
74.60.Jg,74.60.Ge,74.70.Dd 
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Recently, there has been a renewed interest in the study of magnetic superconductors 
following the discoveries of the magnetic superconductor ErNi 2 B 2 C, which, in the 
superconducting phase, has a vortex line-lattice with very unusal properties ||. Ng and 
Varma || proposed that many of the unusual properties of ErNi 2 B 2 C can be understood 
by assuming that the system has an instability towards forming a spontaneous vortex phase 
at low temperatures below Twf ~ 2.3K. The coexistence of weak ferromagnetism and 
superconductivity in the system in the presence of external magnetic field was confirmed in 
neutron diffraction measurement || and in magnetic hysterisis measurement Q recently. 

The studies of Ng and Varma focused at the temperature range T > Twf, when 
long-ranged magnetic ordering is absent. In this paper we shall study numerically the 
H — T phase diagram of a weak magnetic superconductor, where the magnetization from 
the magnetic moment is not strong enough to support a spontaneous vortex phase in absence 
of magnetic field. The behaviour of the system at both T > Twf and T < Twf will be 
studied. In particular, we shall investigate the competition between the spiral phase |§-|8| 
and spontaneous vortex phase. 

We start with the Ginsburg-Landau Gibb's free energy functional G describing a mag- 
netic superconductor system in the presence of external field H, 

+ia|M| 2 + i,3|M| 4 + i 7 2 |VM| 2 - B.M - M], (l) 

where B = Vx A is the total magnetic field, M is magnetization and if) is the superconducting 
order parameter. We assume here that the superconducting component is isotropic in space 
but the magnetic component is restricted to lie only on the x — z plane, as is believed to be 
the case of the ErNi 2 B 2 C compound ]!]. The anisotropy within the x — z plane is neglected 
in this study. The effect of in-plane anisotropy will be discussed at the end of the paper. 
We shall work in the London limit where contributions from vortex cores will be neglected, 
and shall consider only fluctuations in the phase of the superconducting order parameter 
if). In this limit, we obtain an effective GL functional G e ff in terms of the magnetization 
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component and the vorticity (vortex density) field k(x) ||, where 

G eff [M, k]=J d 3 x{~p + i-(£) 2 - ±-Htft + AttM) + \a\M\ 2 + \(3\M\' + \l\ VM| 2 } (2a) 
~ J d 3 xd 3 x'{K -H - 47rM) M (f)(ir 1 )^(£ - x')(k -H- 4ttM) u {x' ), 

with fourier transform of (n _1 ) pi/ (x) given by 

^^^((rri)^-^ 1 )- (2b) 

p = 1^1 2 and Ao are the superfluid density and London penetration depth, respectively, in GL 
theory in the absence of the magnetic component. Notice that only transverse magnetization 
(V x M 7^ 0) couples to vortices. We shall study the spontaneous vortex phase and spiral 
phase of the system in this paper using approximate, variational solutions of G e ff. We first 
consider the case when the external magnetic field H is along z-axis, where the magnetization 
couples directly to H. 

We consider trial solutions of G e ff of form 

M(r, z) = M z + M 1 cos{Qy)z + M 2 sin(Qy)x, (3) 

i 

where r — (x, y). The solution represents an ideal flux-line lattice in z direction with straight 
vortex lines, where Ri = (Xj, Yj) represents the position of vortices in the x — y plane which 
form a regular lattice. The magnetic component may form a spiral state with magnetization 
lying on the x — z plane, where Q = Qy is a spiral wave vector yet to be determined, or form 
a spontaneous vortex phase with M Q ^ 0. Notice that we have also included the possibilities 
of co-existence phases (M , Mi, M 2 ^ 0)in our trial solution. 

Substituting Eq. (|) into Eq. (Q) we obtain after some algebra 



_ = _(M 2 + ^— * + - + M 2 - — ^ + 2M 2 M 2 + ^— i ^ + — i — ^ (4) 



2^ 2 4tt ° 7 8vr 1 + \ 2 Q 2 

1 v B 2 BMi 
8tt f N (1 + \ 2 q 2 N ) (1 + XlQ 2 ) ' 



where V is the volume of the system. The fourier transform of vortex density field n{q, q z ) 
is 

it {0) (q,q z ) = 5{q z )zB5 (2 \q- q N )V, 

where gVs are reciprocal lattice vectors and B = n$ is the average magnetic field trapped 
in the flux-line lattice, n is the density of vortex lines. Notice that M\ couples to the flux-line 
lattice when Q = qj^ for some q^, reflecting the fact that the existence of a flux-line lattice 
induces fluctuations in magnetization which is commensurate with the lattice. 

To simplify our calculation we replace the sum S = ^ 0M g n ^t 2 i \ by the approximate 



value S = BH c i/An |10] which is valid at H ~ H&, where H& is the lower critical field in 
the absence of the magnetic component. Logarithmic correction to S arises at intermediate 
value of H c2 > H > H c i which we shall neglect in the following since we shall be interested 
mainly at qualitative features at regions around H ~ H&. Minimizing the energy with 
respect to Q and assuming that Q ^ q^ we obtain Q~ 2 = Q~ 2 ~ 5^Aq, where Q s is the 
wave vector characterizing spiral instability ||[7| in magnetic superconductors, 5 = ^/AtcXq. 
To examine the importance of magnetization fluctuations commensurate with the flux-line 
lattice, we have also consider Q = q 1 ~ (BIti(k) / H cl )^ Ag 1 in our study, where qi is the 
smallest reciprocal lattice vector of the flux-line lattice. The corresponding free energies are 
minimized with respect to B, M D , Mi, M 2 at various temperatures and external magnetic 
fields H to determine the phase diagram. We have chosen parameters such that temperatures 
are measured relative to superconducting transition temperature T c , which is set to be 
T c = 10. We have chosen also k = Ao/£ = 5, and a = a (T/t m — 1), where a Q = 10 
and (3 = 20. t m is chosen such that Tm = 3 (or 0.3T C ), where Tm is the ferromagnetic 
transition temperature in the absence of superconducting component. Notice that in the 
presence of the superconducting component, the actual magnetic transition temperature 
Twf is below Tm because the Meissner effect forbids appearance of uniform magnetization 
in a superconductor |5|,[7J . Moreover, the zero field magnetization M sat at T = in absence 
of superconducting component is chosen to be M sat = H c ijAm, i.e. the system is marginal 



in supporting a spontaneous vortex phase. 

Minimizing the energy of the system we find that magnetic fluctuations commensurate 
with flux-line lattice is always unimportant compared with the spiral instability and the 
system is in a spiral state with Mi = M 2 and M Q = at zero magnetic field at temperatures 
T < T WF ~ 2.74. The system goes through a second order phase transition to a non- 
magnetic superconducting state at T = Twf- At T < Twf, the system stays at the spiral 
state when external magnetic field is applied, until H — > H cs < H ci , where the system 
undergoes a first order phase transition into a magnetic-field assisted spontaneous vortex 
phase with finite M a , and with Mi = and small M 2 . The presence of M 2 ^ is a residue 
effect of spiral instability in the spontaneous vortex state. We find numerically that M 2 
reduces rapidly to zero as external field H increases further. The H — T phase diagram is 
shown in Fig.l, where (VP), (SP), (MP) denotes the magnetic field-assisted spontaneous 
vortex phase, spiral phase and Meissner phase, respectively. H* x is the zero temperature 
lower critical field in the absence of the magnetic component. The dotted and dashed lines 
give H c2 and H cl of the superconducting component, respectively, in the absence of the 
magnetic component. The solid line gives H cs and represents a first order transition line 
from the spiral phase to the vortex phase at T < T WF , and represents the usual transition 
from Meissner phase to vortex phase at T > Twf- The dot-dashed line denotes the second 
order transition from Meissner phase to spiral phase at H < H cs . 
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FIG. 1. H-T phase diagram of the magnetic superconductor 



To make comparisons with experiments we show in Fig. 2 the total magnetization M 
(B — H)/An versus external magnetic field HjH.* x at three different temperatures T 
0.0,2.0,4.0, 
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FIG. 2. Magnetization versus magnetic curve for at three different temperatures T=0.0, 2.0, 4.0. 

The qualitative behaviour of the M — H curve is similar to those observed experimentally 
in the ErNi 2 B 2 C compound at low temperature ||. Notice a sudden jump in magnetization 
at H = H cs appears when T < Twf, signaling the first-order transition to magnetic field 
assisted spontaneous vortex phase. This sudden jump is absent at T > Twf, justifying 
the identification of Twf as the spiral state transition temperature. The sudden jump in 
magnetization will probably be smeared out by disorder in realistic samples, and H cs can 
be considered as an effective lower critical field observed in magnetic superconductors at 
T < Twf- Notice that in previous analysis of Ng and Varma ||, the lower critical field H c \ 
for magnetic superconductors was computed from single vortex line energy and it was found 
that the value of H c i is not much affected by the magnetic component in the spiral state. 
However, we find here that H cs is smaller than H c \ at T < Twf-, indicating that an effective 
reduction in H c i can occur because of the first order transition to the magnetic field assisted 
spontaneous vortex phase. As M sat increases further, we find numerically that H cs starts 
to decrease as temperature decreases until when AnM sat ~ 4iJ* x where H cs — > at T — > 0, 
indicating a zero temperature phase transition of the system to the spontaneous vortex 
phase at zero magnetic field. As M sat increases further, the region where the spontaneous 



vortex state is stable at zero field expands further and the spiral state becomes stable only 
at a finite temperature range < T sp < T < T WF , in agreement with results from previous 
studies [7p. 

Recently it was observed in magnetization studies of the ErJS 'i^B-xC compound that 
substantial hystersis developed in the M — H loop at temperature T < T WF ||, and the 
result was interpreted as coming from a ~ three fold increase of the pinning force of the 
flux line lattice at T < Twf- Here we shall consider the effect of spiral phase in magnetic 
hystersis. We note that the spiral phase and spontaneous vortex phase are locally stable 
phases at T < Twf over a wide range of magnetic fields. The spiral phase is (locally) 
stable at H < H[, and the spontaneous vortex state is (locally) stable at H > H l 2 , where 
H\(H l 2 ) depends on the parameters in the GL functional, and is a decreasing (increasing) 
function of temperature, with H[ > H l 2 in general. We propose that the huge magnetic 
hystersis observed at T < Twf is a result of the system being trapped in the two meta-stable 
states when magnetic field increases (decreases) at temperature T < T WF . When magnetic 
field increases from zero, the system is trapped in the spiral state until H > H±, where 
the spiral state becomes locally unstable and the system goes to the spontaneous vortex 
state. Similarly the system is trapped in the spontaneous vortex state when magnetic field 
decreases, until H < H l 2) where the system goes back to the spiral state. As a result, huge 
magnetic hystersis can be observed at magnetic field range H l 2 < H < H\. 
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FIG. 3. Magnetic hystersis between spiral and (magnetic-field assisted) spontaneous vortex 
phase at three different temperatures T=0.0, 1.1, 2.2. 

In figure 3 we show the hystersis magnetization AM = \M up — M<} own \ as function of 
external magnetic field H/H^ at three different temperatures T = 0.0, 1.1, 2.2, assuming the 
above picture. The narrowing of the hystersis region and reduction in AM as temperature 
increases is in qualitative agreement with what is observed experimentally [f|]. Notice that 
first order jumps of AM is found at H = H[ and H = H\ in our theory. These jumps should 
be smeared our by defects in real samples. Notice also that additional features are observed 
experimentally at H ~ Hf at low temperatures |4[], suggesting that there exist additional 
structures in the spiral-spontaneous-vortex state phase diagram not included in our simple 
variational solution. 

Next we consider the case when the applied magnetic field is almost perpendicular to 
the magnetic plane and makes a small angle 9 with the ^-axis. In this case, the coupling 
between magnetization and external field is weak. However, it was observed experimentally 
that the flux-line lattice seen to rotate away from the applied field direction at temper- 
ature T — > T WF [Q. This observation can be explained qualitatively by development of 
in-plane magnetization as T — > Twf Hi- Here we shall study numerically this tilting effect 
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at temperatures both above and below Twf- 

The trial solutions of G e // we consider is similar to Eq.(3), except that /t(r, z) is in Eq.(3) 
is replaced by 

k(t,z) = $ J2(cos(6 t )y + sin(9 t )z)5 {2 \r - % 

i 

where 9 t is the angle the flux-line lattice made with the y-axis. Correspondingly, the term 
-B(H + 4vrM )/4vr in Eq.(4) is replaced by 



- B ^ H \ 4?rMo) -> -^B {Hcos{9 - 9 t ) + 4KM o sin{0 t )) . 

47X 4:71 

Minimizing the energy of the system we obtain tan(9 t ) = (Hsin(0) + 4ixM ) / Hcos{6). 
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FIG. 4. Angle of the flux4ine lattice with the y-axia as a function of magnetic field at three 
different temperatures T = 0.0,1.4,2.4. The angle between applied magnetic field and y-axia is 
1.8°. 



The resulting angle 9 t is shown in figure 4 as a function of external magnetic field H/ H* x 
at three different temperatures T = 0.0, 1.4, 2.8 with 6 = 0.0br(1.8°). At T = 0, 9 t remains 
at 1.8° at low field, indicating that the system is in the spiral state. Flux-line lattice with 
9 t = 9 (M Q = 0) begins to form at H ~ H cl , because the coupling of the magnetization to 
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the applied field is very weak at small 9. A first order phase transition to the (magnetic-field 
assisted) spontaneous vortex phase occurs at H ~ 1.6if* l5 where in-plane magnetization M 
becomes nonzero (Mi = 0, M2 small) and 8 t jumps to a value ~ 25°. As magnetic field 
increases further, 9 t decreases because the in-plane magnetization is already near saturation 
and the rate of increase in M is slower than the rate of increase in H . Similar result is 
also found at T = 1.4, except that the first order transition occurs at a higher field and 9t 
is smaller. We find numerically that there exists also a small region between the spiral and 
spontaneous vortex phase where the two states are co-existing in the system with a small 
value of M but large Mi ~ M2, with 9 t ~ 4°. This state is stable only at intermediate 
temperatures below T WF and at a narrow range of magnetic field H > H cl . At temperatures 
above but close to Twf, we find a continous increase of 9 t as function of H at H > H c \, 
in agreement with previous experimental observation ||. When 6 = 0, direct first order 
transitions to spontaneously-tilted flux-line lattice phase is also found at T < Twf, when 
magnetic field is large enough. 

Summarizing, using simple variational solutions of the GL free energy functional, we 
have investigated numerically the magnetic behaviours of a magnetic superconductor where 
the magnetic component is weak and marginal in supporting a spontaneous vortex phase. 
We have studied both the cases when magnetic field is parallel to the magnetic plane, 
and when it is almost perpendicular to the magnetic plane. Our results are in general 
agreement with experimental observations. Furthermore, we have made several theoretical 
predictions in this paper, including the appearance of spiral phase at zero field at T < T WF , 
which can be tested by neutron scattering experiment. Notice that we have neglected in- 
plane magnetic anisotropy in our study. In the case of ETN12B2C where magnetization are 
strongly confined to point at either x or z directions, a smooth spiral phase is impossible 
and will be replaced by sharp domain structures H that can be disordered easily. In this 
case the sharp elastic neutron peaks at spiral wavevector Q s will be replaced by a broadened 
peak centered at wavevector q = with width of order Q s . We predict also a first order 
transition to spontaneous vortex phase when external field is strong enough, both when the 
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field direction is in the magnetic plane, and when it is almost perpendicular to the plane. 
In general a stronger magnetic field is needed for this transition to occur when the angle 
between the applied field and the magnetic plane is larger. When the magnetic field is almost 
perpendicular to the plane, the first order transition to spontaneous vortex phase appears as 
a spontaneous tilting transition of the flux-line lattice, which can be tested experimentally. 

T.K. Ng acknowlegdes many helpful discussions with CM. Varma. This work is sup- 
ported by Hong Kong UGC through Grant No. HKUST6124/98P. 
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